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Abstract 

The state space of a polymer molecule is anal- 
ysed. We show how the size of the state 
space can be reduced on the basis of sym- 
metry. In the reduced state space, the proba- 
bility of a new state (termed below as class) is 
equal to the number of old states represented 
by the new state multiplied by the proba- 
bility of each old state. As an application, 
the electrophoretic motion of the molecule in 
gel is considered. We discuss the influence 
of the gel medium and of external field on 
the molecule states, with absorbing states of 
hooked molecules playing a major role. We 
show that in the case of strong fields both the 
velocity and the diffusion coefficient decrease 
with field. Finally, we evaluate the time of 
relaxation to and from the absorbing states. 
This is done with a continuous version of the 
exact enumeration method for weighted net- 



works. 



1 Introduction 

Computational tasks in biology, sociology, 
physics and many other areas meet in most 
cases limitations of computer memory and 
time, which limits the size of the modelled 
systems. It is worthwhile to look for methods 
how to reduce the size of an analysed system 
without loss of information. Recently we 
have shown [H [2] that when some symmetries 
are present, the size of the state space can 
be significantly reduced. In our previous 
papers we applied the new method to spin 
lattices [I] and traffic systems [2]. In both 
cases the same formalism was used. The 
idea is to group similar states on the basis 
of their symmetry. The procedure preserves 
information on the static probabilities of all 



states. 

The same idea is now applied to the state 
space formed by a set of possible confor- 
mations of a circular polymer molecule. 
A standard model description of this kind 
of molecules is the repton model [3l H]. 
The main trick of this model is to divide 
the whole molecule into parts, so-called 
reptons. We are interested in all possible 
conformations of a molecule, assuming that 
reptons are placed at nodes of a regular, 
square lattice (another basic assumption of 
the repton model). This means that the 
distance between two neighbouring reptons 
can be equal to or 1 in units of lattice 
constant. In other words, two neighbouring 
reptons can be placed in the same position 
or in an assumed, constant distance. A 
symmetry of the space will also be visible in 
molecule conformations, which allows for a 
reduction of the size of the state space of the 
system. 

Modelling of conformations of polymer 
molecules and transitions between them 
finds applications, e.g. in electrophoretic 
transport in gel. This process is commonly 
used in biology for separating molecules, 
such as DNA or peptides, according to their 
size. The process is based on the motion of 
charged polymer molecules in an external 
electric field in a porous medium. In the 
theory of electrophoresis one can distinguish 
different regimes of molecule migration 
through the media, which depend on the 
molecule length, the gel density and the 
electric field strength. In the case of weak 



fields the process of molecule displacement 
is a result of thermal diffusion and motion 
caused by the applied electric field. In the 
case of strong fields, the picture changes 
significantly as thermal shifts become more 
rare. It may lead molecules to hook on gel 
fibres, which reduces their speed or even 
immobilize them |5l |6] . 

The paper is organised as follows: First, 
the state space of the molecule conformations 
is presented. We show how the symmetry 
of the circular molecule allows the reduction 
of the size of the system. Then we present 
classes of states for the original and the re- 
duced system. In the following section we 
show how the state space is changed because 
of the presence of the porous medium and an 
external electric field. Section [3] is devoted 
to the discussion of the probability of states 
and the probabilities of particular transitions 
between states in an external field. The fol- 
lowing section presents the results of our cal- 
culations of the velocity of molecules, the dif- 
fusion coefficient and the relaxation time of 
the system as dependent on the applied field. 
In the case of the relaxation time, the reduc- 
tion of the state space leads to approximate 
solutions. The final section is devoted to dis- 
cussion. 



2 Analysed system 

2.1 State space 

The state space is formed from the set of all 
possible conformations of a circular polymer 



molecule in the repton model. We also as- 
sume that the reptons are indistinguishable. 
This means in particular that a shift of the 
whole molecule along its contour does not 
lead to a different state. We note that the 
lattice cells are distinguishable. However, 
a translation of the whole molecule in the 
lattice does not generate a new state. With 
these assumptions, we find all possible con- 
formations of a molecule of a given length. 
The length of the molecule is expressed by 
its number of reptons A^. 




2.2 Network of states 

A transition between two molecule conforma- 
tions is possible if it is connected with the dis- 
placement of one repton from one lattice node 
to another one. One can think about states 
as about nodes of some network of states, 
with edges determined by possible one-repton 
moves between the states. In this terminol- 
ogy the number of connections of a given 
state to other states is equivalent to the node 
degree. The network obtained for a circular 
molecule of five reptons is presented in FigjT} 
In Fig|2] all possible transitions from one 
example state are presented for a molecule of 
5 reptons. All of them are connected with the 
change of the position of one repton, and the 
distance between two neighbouring reptons is 
equal to the lattice constant or to 0. 



Figure 1: Network of states for 5-repton long 
molecule. 




Figure 2: Examples of possible one-repton 
transitions between states. 



2.3 Classes 

As points in the space are distinguishable, 
molecule conformations which are mirror re- 



flections are seen as different states. However, 
in fact tlie distinction between tfiose states 
is illusory. Existing symmetry of tlie space 
causes tliat all states which can be obtained 
one from another as a result of reflection can 
be treated as one class. As a consequence, 
the effective size of the network of states can 
be reduced. 

Formally the procedure of reduction of the 
size of the state space of the system is per- 
formed as follows: The quantity which char- 
acterises each state is the number of transi- 
tions between this state and other states of a 
system. The number of possible transitions 
from a given state and the set of states it can 
be transformed to determine a class to which 
the state is classified. Each conformation of 
the molecule is then represented by the class 
it is classified to (see Figjs]). 
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Figure 3: Examples of possible one-repton 
transitions between classes. 



The inspection of all possible transitions 
leads to the division of the whole set of states 
into classes. The size of a state space and the 
number of classes depend on the length of the 
molecule. Obtained results are presented in 
TabjT] As it can be seen in these data, the 
number of possible states grows significantly 
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NUMBER OF STATES 


NUMBER OF CLASSES 


5 


35 


9 


6 


123 


31 


7 


359 


62 


8 


1363 


241 


9 


4 891 


702 



Table 1: Number of states and number of 
classes for molecules of different lengths A^ 
(expressed via the number of reptons) 
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CLASSES GIVEN CLASS 


CAN BE TRANSFORMED 


TO 


DEGREE 


#STATES 


Al 


Bl 


Bl CI El 
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A2 


E2 


E2 E2 E2 
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Bl 


Al 


B2 CI Dl 


Fl 






5 
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B2 


Bl 


Bl B2 B2 


Fl 






5 


4 


CI 


Al 


Bl Bl CI 


Dl 


E2 




6 


4 


Dl 


Bl 


Bl CI Dl 


E2 


Fl Fl 




7 


4 


El 


Al 


Al E2 E2 


Fl 


Fl Fl 


Fl 


8 


2 


E2 


A2 


CI Dl El 


E2 


E2 Fl 


Fl 


8 


4 


Fl 


Bl 


Bl B2 Dl 


Dl 


El El 


E2 E2 


9 


4 



Table 2: Classes identified for the 5-repton 
long molecule 

with the number of reptons, so in the follow- 
ing part of the paper we will work with a 
molecule of 5 reptons. 

A full classification of the states of a cir- 
cular molecule of 5 reptons is presented in 
Tab|2] The symbol in the first column indi- 
cates a class. In the second column classes 
of states are listed which are attainable from 
the states classified to the class specified in 
the first column. The third column covers in- 
formation about the number of possible tran- 
sitions, and in the last column the number of 
states which are classified to a given class is 
given. 

According to the Tabjlj for 35 possible con- 
formations of the 5-repton long molecule, 9 
classes allow for the correct description of the 
system. The probability of a particular class 
is equal to the sum of probabilities of states 
classified to this class. Here again, the ob- 
tained system can be described in terms of 



networks: nodes are identified with classes, 
and ties - with connections between states 
belonging to a given class. Because of that, 
the obtained network is weighted. The net- 
work of classes is presented in Figji} 




Figure 4: Network of classes of 5-repton long 
molecule. 



2.4 Molecule in porous medium 

The first consequence of the presence of the 
medium is that not all theoretically possi- 
ble molecule conformations are possible. An 
example is presented in FigjS} where one of 
the possible conformations of the 5-repton 
molecule is shown. The conformation pre- 
sented in panel a (and three similar confor- 
mations resulting from the space symmetry) , 
is allowed. However, if we analyse molecules 
moving in the porous medium, such confor- 
mations cannot be obtained. The molecule 
encircles the fibre in the middle of panel h of 
FigjSJ As we assume that both ends of each 
fibre are connected with the medium, this 



situation cannot happen. Because of that, 
some configurations have to be removed from 
the state space. The number of conforma- 
tions which should be removed depends on 
the length of the molecule. In the case of 
five reptons, the four earlier mentioned con- 
formations cannot occur. Because of that, 
the state space contains now 31 states. What 
is also changed is the set of possible transi- 
tions between states, as the fibres of a gel act 
for migrating molecules as obstacles. From 
the originally specified set, we have to re- 
move the transitions which pass across a fi- 
bre. An obvious consequence is the change of 
the number of classes. The result is presented 
in Tabm 



b) 



n 



X 



Figure 5: a) One of the possible con- 
formations of the 5-repton molecule. b) 
Some conformations becomes impossible if 
the molecules are embedded in a gel medium 
(the polymer is linked to a fibre with cross- 
section denoted by the solid bullet in the mid- 
dle). 



The number of classes which initially equal 
to 9 is reduced to 6. The minimal degree is 
now equal to 2, and the maximal one equals 
6. It results from the fact that now not 
only are some conformations not allowed but 



CLASSES GIVEN CLASS CAN BE TRANSFORMED TO DEGREE #STATES 



Al 


CI 


Dl 
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12 


Bl 


Bl 


Dl 


D2 
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CI 


Al 


Al 


D2 


D2 
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C2 


D2 


D2 


D2 


D2 






4 
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Dl 


Al 


Al 


Al 


Bl 


Dl 


D2 
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D2 


Bl 


CI 


CI 


CI 


C2 


Dl 
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4 



# LINKS 


#CLASSES 1 


#CLASSES < 


2 


7 


4 


3 


3 


2 


4 


5 


4 


6 


6 


4 



#STATES 



12 
4 

7 



Table 3: Classes identified for the reduced 
state space of a circular molecule of 5 reptons 
in gel medium. 

some repton displacements are also prohib- 
ited. Looking at Fig|3]we can see that tran- 
sitions from the conformation in the middle, 
classified originally to the class El, to the 
states classified to class Fl become impossi- 
ble in the presence of gel fibres. As a conse- 
quence, the degrees of states and the number 
of classes decrease. One class from the orig- 
inal set of classes disappears as it contains 
states which are now impossible to obtain. 
Classes Al and 51 of the original system are 
now merged, and the degree of states belong- 
ing to a new class is equal to 2. Also original 
classes El and Fl are merged. Other classes 
remain unchanged with respect to states they 
include; what does change is the degree of 
states. The only exception is one-element 
class which consists of the state in which all 
reptons are placed in the same lattice cell. 
This class is identified in both cases. The 
number of possible transitions from this state 
in both cases is equal to 4. 

2.5 Molecule in porous medium 
and electric field 

The presence of an external electric field 
causes sjmametry of the system to be partially 



Table 4: Classes identified for field directions 
parallel (||) and diagonal (<) to the lattice 
axis. 

classes becomes higher than in the case with- 
out an electric field. In this case, the network 
of states is directed, with weights of transi- 
tions between given two states determined by 
the actual value and direction of the applied 
field. The class of states is now specified by 
their out-degree. Although, the number of 
links going in and out of each state is the 
same, their weights are different. 
The number of classes depends on the direc- 
tion of the applied field with respect to the 
square lattice. The results are presented in 
Tab|4j If the direction of the field is par- 
allel to the lattice axes, 21 classes for the 5- 
repton long molecule are identified. A slightly 
smaller number of classes - 14 - is identified 
for directions parallel to the diagonal of the 
lattice cells. 

3 Probability distribution 
dynamics 

3.1 Motion in field 

The possibility of the molecule's motion de- 
pends on the electric to thermal energy ratio 
and the direction of the applied electric field. 
Those two agents determine the probabilities 
of all possible shifts of the molecule in a given 



lowered. As a consequence, the number of conformation. In the case of motion in the 



strong electric field, shifts which cause an in- 
crease of energy become practically impossi- 
ble. The probabilities of those transitions are 
very close to zero. In this case, it is possible 
that for a given field direction some confor- 
mations may behave as traps. This situation 
occurs for the states for which all possible 
transitions are connected with a displacement 
of the repton in the direction opposite to the 
direction of the applied external electric field. 
A change of the molecule's position in such a 
conformation requires a very long time. We 
observe also that some states which are not 
traps themselves, quickly achieve a trap con- 
formation as a result of the allowed transi- 
tions. However, for a majority of states a 
series of transitions leading to the molecule 
displacement in the space are possible. For 
moderate and weak fields, trap states are not 
observed, in the sense that probabilities of all 
potentially possible transitions are non-zero. 
In all cases, the probabilities of transitions 
are calculated in accordance with the for- 
mula: 



w{e) 



exp{- 



-ecosa) 



(1) 



3.2 Stationary probability dis- 
tribution of the molecule 
conformat ions 

In the case of the absence of any external 
force, the probabilities of all states are the 
same. Switching an external electric field on 
causes the differences in probabilities which 
depend on the field strength. For each case a 
transition matrix can be constructed. How- 
ever, the problem may occur in the case 
of strong fields. The existence of absorb- 
ing (trap) states causes that corresponding 
columns in the transition matrix contain ze- 
ros. In this case, the probabilities of the 
states depend on the initial state. Taking 
this into account, we use the set of Master 
equations [7]: 

dMi{t) 
dt 

where for a given state i summing up goes 
through its neighbours Si with probabilities 
of transitions w calculated for a given value 
and sense of the external electric field (in ac- 
cordance with Eqjl]). The stationary distri- 
bution is found as the asymptotic one for long 
times, with the probabilities of all allowed 
states initially equal. 



> ]M,{t)w,^,- 


-^M,{t)wi^j 


jes. 


jes. 




(2) 



where: e = AE/kT - electric to thermal 
energy ratio, a - angle between the currently 
considered electric field vector and the 
direction in which a given repton changes its 
position. 



3.3 Exact enumeration and 
its continuous version for 
weighted networks 

Trap conformations may not be observed for 
all molecule sizes and not for all field direc- 
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tions. In the case of a 5-repton molecule, 
trap conformations are observed if the ex- 
ternal field is directed along the diagonal of 
the lattice cells. For each of these four possi- 
ble field directions, two trap conformations 
are identified. In this case one can evalu- 
ate the time after which, starting from any 
non-trap conformation, the molecule will be 
in the trap state. To find this time we use 
an approach similar to the exact enumera- 
tion method (EE) 0. 

The advantage of this method is that it allows 
for taking into account the whole set of initial 
conditions at the same time. In EE, a record 
is constructed of length equal to the number 
of possible states of the system. In the initial 
state, values for all non-trap states are set to 
one, and for the trap states to zero. Each el- 
ement of the record keeps the information of 
the probabihty of a given state. At each step 
of EE procedure for all non-trap elements, a 
new value is calculated, while elements for 
trap states are for the whole time kept un- 
changed (with value equal to 0). The value 
currently assigned to each state is equally di- 
vided into all its neighbours. So a new value 
of each state is calculated as a sum of the 
fractions of probabilities obtained from states 
which can be transformed to this state. An 
important assumption in [8] is that all tran- 
sitions are equally probable. 
In our case this assumption cannot be main- 
tained. Probabilities of transitions to differ- 
ent states are of different values and they 
change with the external field direction and 
strength. To take this into account, here we 
propose a slightly different approach. Sim- 
ilarly to the original method, we construct 



a record of the size equal to the number of 
states. The initial values are set to 1 for 
non-trap elements, and to for trap ele- 
ments, with the latter values kept unchanged 
for the whole time. The difference relates 
to the method of calculating the new value 
for the non-trap states, which allows taking 
into account different probabilities of partic- 
ular transitions. The idea is to use Master 
equations (Eq|2]). Such an approach reflects 
the different time scales of different transi- 
tions. The method proposed is a continuous 
version of the exact enumeration procedure 
for weighted networks. Below it will be used 
to evaluate the relaxation time to and from 
the absorbing states. 

4 Results 

Both the velocity and the diffusion coefficient 
are found from the stationary probability dis- 
tribution of classes. There, the probability of 
each class is equal to the number of states rep- 
resented by the class multiplied by the prob- 
ability of each of these states. On the con- 
trary, the results on the relaxation time are 
approximated. Below they are presented for 
the original and the reduced state space. 

4.1 Velocity of the molecules in 
field 

For molecules migrating in the porous 
medium one can estimate the velocity of dis- 
placement as a function of e. To obtain this 
quantity we calculate for each value of e the 
probability of each molecule conformation in 



the equilibrium state Pi{e). For this purpose 
we use the set of Eqsj2] and we find numeri- 
cally the solution in asymptotically long time, 
with an initial state where each conformation 
is equally probable. Then we calculate the 
field dependence of the velocity, using the for- 
mula: 



N 



v{e) = ^Pi{e) X {wi^j{e) x k) 



j&s^ 



(3) 



where: Pi{s) - probability of a given state at 
a given value of e, Wi^jie) - probability of a 
transition from state i to state j at a given 
field e, k G (—1,0, 1) - direction of a repton 
displacement in relation to a given external 
field direction. 

The obtained result depends on the con- 
sidered direction of the external electric field. 
The reason is that different behaviour is ex- 
pected depending on the existence or lack 
of the absorbing trap conformations. Their 
presence will cause a decrease of the molecule 
velocity in the case of a strong electric field, 
as their probability in the equilibrium state 
then becomes dominant. Results for the 5- 
repton long molecule are presented in Figj6] 

The line with squares presents results ob- 
tained for the case free from the trap states, 
i.e. when the field is parallel to the lattice 
axes. In this case the velocity grows linearly 
with the growth of e for values of e ^ 2. Then 
the curve begins to saturate and reaches con- 
stant value for e ^ 5. 

If the molecules can get trap conformations, 
the results are qualitatively different. For 
€ ^ 1 the curve (line with dots in FigJGJ) 
behaves approximately in the same way, as 




Figure 6: Velocity as a function of e for the 5- 
repton long molecule. The line with squares 
presents results for the displacement direc- 
tion parallel to the lattice axes. The line with 
dots presents results for the displacement di- 
rection at an angle 45° to the lattice axes. 
Here e is unitless, temperature is set to 1, 
and the velocity is expressed in the number 
of cells per computer time step. 
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for the field direction parallel to the lattice 
axes. In this range of e, an increase of the 
velocity with field is observed. The influ- 
ence of the trap states becomes visible for 
e ^ 2, when a maximum of the velocity is 
found. For e G (1, 10) the velocity decreases, 
as exp(— ae). Accurate calculation of the a 
parameter is difficult due to the shape of the 
curve but it can be estimated at ~ 0.5. In [9], 
the same character of the field dependence of 
the velocity was found for linear molecules, 
with the coefficient a = 3/2. For e ^ 10 the 
calculated velocity achieves zero with numer- 
ical accuracy. 

4.2 Diffusion coefficient in 
the presence of absorbing 
states 

Apart from the velocity, the diffusion coeffi- 
cient as a function of e can also be estimated 
from the probabilities Pi{e) in an equilibrium 
state. We use the formula: 



D{e) oc 



iatrap 



E P^^ 
idnot trap 



ficient D is very small and almost constant for 
e < 0.1, then increases linearly for e G [0.1,2] 
and achieves maximum for e f« 3. In the 
range e G (3, 10) the diffusion coefficient de- 
creases to 0. Similarly to the velocity depen- 
dence, the value of the decrease parameter 
can only be estimated and is ~ 0.43. The 
shape of the curve is consistent with the ve- 
locity relation in the presence of trap states. 




(4) 



Figure 7: Diffusion coefficient Z^ as a func- 
tion of e for the 5-repton long molecule in 
the presence of absorbing states. Here e is 
unitless, temperature is set to 1, and D is in 
lattice units squared per computer time step. 



In other words, the diffusion coefficient is 
calculated as a product of the probability 
that the system is in the trap state (the first 
factor) and the probability that the molecule 
migrates in a media (the second factor) . This 
is an approximation of the binary distribution 
with two states, where molecules move with 
constant velocity in one state and they are 
trapped in another state. 

The obtained results are presented in FigjT] 
in semi-logarithmic scale. The diffusion coef- 



4.3 Relaxation time to and 
from the absorbing states 

To estimate the relaxation time of the 
system we calculate the average value of M 
of the whole set of non-trap states at each 
step of evolution. The method uses Master 
equations and the exact enumeration tech- 
nique, as described in the last part of Section 
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3.3. As the result, we obtain the curve 
which decreases as exp{—t/T). Parameter r 
estimates the time needed to achieve trap 
states, against e. The same formahsm can 
be used to find the time needed for an escape 
from the trap conformations, provided that 
this process is possible in a finite time. As 
it was discussed earlier, this is possible if a 
shift caused by the thermal energy can occur 
against the field; this is the case in moderate 
and weak electric fields. In the case of a 
strong field, the molecule which achieves a 
trap conformation will remain in this state 
practically forever. 




10 100 



Figure 8: Time of enrichment the trap state 
Ttrap (left hand side vertical axis) and time 
of escape from the trap Tnon trap (right hand 
side vertical axis) as a function of e. The 
empty symbols refer to the states network, 
while the full one to the classes network. Here 
e is unitless, temperature is set to 1, and r is 
in computer time steps. 



The obtained results for a 5-repton long 
molecule are presented in FigJs] in semi- 



logarithmic plot. The lines with dots show 
time Ttrap (left vertical axis) needed to achieve 
trap state (class of traps) as a function of e. 
The empty symbols refer to the network of 
states, and the full ones to the network of 
classes. In a very low (below ^ 0.1) and in 
a high (above ~ 5) range of e, Ttrap does not 
change with the change of the field strength. 
In the first case the time of trapping is long, 
which is connected to possible thermal move- 
ments of reptons. In the case of a strong 
field it is clear that molecules will be quickly 
stopped in trap configurations. In between 
the two constant ranges the curve has a steep 
slope. An increase of the external electric 
field makes thermal shifts impossible, which 
shortens the time of trapping. 
The line with squares shows time Tnon trap 
(right vertical axis) needed to escape from 
trap states (class of traps) as a function of e. 
For values of e < 1 the escape from the trap 
conformation occurs in a short time. When 
the field increases, the escape time begins to 
increase quickly and without an upper limit. 
For both curves, a qualitative agreement of 
results achieved for the network of states and 
the network of classes is obtained. Shorter 
relaxation times obtained in the case of the 
network of classes are the consequence of re- 
ducing the system size. 



5 Discussion 

We have demonstrated that the formalism 
of the reduction of the size of the state 
space of systems based on system symmetry 
can be applied in the case of state space of 
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molecule conformations. The rate of the size 
reduction depends on the molecule length, 
and it is also different in the presence a of gel 
medium and the external electric field. In 
the last case the reduction of the system size 
is not as significant as in the absence of a 
field, because of the field-induced symmetry 
breaking. A qualitative agreement of the 
results of the relaxation time for the orig- 
inal network of states and for the reduced 
system of classes is obtained. It shows that 
the system properties are preserved in its 
compressed representation. 

For the obtained states we show how 
probabilities of states can be calculated. The 
presented approach can be applied for any 
weighted system, even if absorbing states 
are observed in the system. We also show 
how the exact enumeration method can be 
generalized to be used in the case of weighted 
networks. 

Although the results are presented for 5- 
repton long circular molecules, the proposed 
approach can also be applied for different 
lengths and shapes of molecules. Also in 
this case, some symmetries will be observed. 
The effects discussed above are expected 
to be present also for molecules of other 
shapes and lengths. In particular, in strong 
fields trap states are generic. Their absence 
in our case of a field parallel to the lattice 
axes seems to be an artifact of the lattice 
itself. The result that velocity decreases 
with the electric field strength seems to 
be paradoxical. However, a similar result 
was obtained in an experiment with DNA 
knots [To] and in another theoretical model 



for the linear polymers migrating in strong 
electric fields [9]. Our results obtained for 
the reduced state space are consistent with 
literature. 

We conclude that the method used to 
reduce the size of state space is general and 
useful numerically. 

Acknowledgement: The author is grate- 
ful to Krzysztof Kulakowski for critical read- 
ing of the manuscript and helpful discussions. 
The research is partially supported within the 
FP7 project SOCIONICAL, No. 231288. 

References 

[1] M.J.Krawczyk, "Symmetry induced 
compression of discrete phase space", 
Physica A 390, 2181 (2011) 

[2] M.J.Krawczyk, "Topology of space of 
periodic ground states in antiferromag- 
netic Ising and Potts models in selected 
spatial structures". Physics Letters A 
374, 2510 (2010) 

[3] P.-G. de Gennes, "Scaling concepts in 
polymer physics", Cornell University 
Press, 1996 

[4] M.E. J.Newman, G.T.Barkema, "Monte 
Carlo methods in Statistical Physics", 
Clarendon Press, Oxford 1999 

[5] J.L.Viovy, "Electrophoresis of DNA and 
other polyelectrolytes: Physical mecha- 
nisms". Review of Modern Physics 72, 
813 (2000) 



12 



[6] G.W.Slater et al., "Theory of DNA elec- 
trophoresis" , Electrophoresis 23, 3791 
(2002) 

[7] N.G.van Kampen, "Stochastic processes 
in physics and chemistry", Elsevier Sci- 
ence Publisher B.V. 1981 

[8] S. Havlin et al., "Exact enumeration of 
random walks with traps". Journal of 
Physics A 17, L347 (1984) 

[9] A.B.Kolomeisky, A.Drzewinski, "Poly- 
mer dynamics in repton model at large 
fields" , Journal of Chemical Physics 120, 
7784 (2004) 

[10] C.Weber et al., "Gel Electrophoresis of 
DNA Knots in weak and strong elec- 
tric fields" , Biophysical Journal 90, 3100 
(2006) 



13 



